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Abstract. 

O , We prove that certain quotients of entire functions are characteristic functions. Under some conditions, 

the probability measure corresponding to a characteristic function of that type has a density which can be 
expressed as a generalized Dirichlet series, which in turn is an infinite linear combination of exponential or 
_ . Laplace densities. These results are applied to several examples. 
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1 Introduction 



There are many cases in the literature where a characteristic function ip(t) of a probability measure can be 
written as (p(t) = l/g(it) for t € R, where g(z) is an entire function of the complex variable z. Two important 
examples are the square of a Kolmogorov law and the Levy area. Specifically, the characteristic function 
of the square of a Kolmogorov law is given by ipx(t) = \/2it/ sin(\/2£i), and thus we are setting gi(z) = 
sin(y / 2z)/v / 2z. This example was studied by Dugue 0, who determined that its distribution function is the 
Jacobi theta function 

oo 

= 4 (x) = 1 + 2 Y^{-l) k e- n2k2x/2 , x>0. 
fe=l 

The Levy area was introduced by Levy lPl4ll (with a slightly different parametrization) as the random variable 
given by 

W 1 (t)dW 2 (t) - [ W 2 (t)dW 1 (t), 



o Jo 



X ■ 

where {Wi(i), t > 0} and {W2(t), t > 0} are two independent standard Brownian motions. Levy |[T4ll 
deduced that its characteristic function is ip%{t) = sech(t). In this case we take g 2 {z) = cos(z). The density 
was also computed by Levy Ifl4ll and is 

/ 2 (x) = -sech(-x) =— r-r = V -1 w e ' W , x^0. 

fc=i 

What these two examples share is that regardless of the behaviour of the entire functions g\(z) and g 2 (z), 
their inverses are characteristic functions. In fact, it is well known that for an entire function g(z) of order 
p < 2 which has only real zeros and g(0) = 1, the inverse tp(t) = l/g{it) is a characteristic function (see 
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Lukacs Ifl5l pp. 88 and 212] for an equivalent formulation). This result follows from Hadamard factorization 
Theorem (see, for example, Levin |[T3l p. 26]) which states that an entire function of order p < 2 can be written 

as 



where {a n , n > 1} are the zeros of g{z) (c must be real so that the characteristic function assertion is true). 
Note that such factorization of g(z) induces a factorization of <p(t), and each factor is of the type 



which is the characteristic function of a translated positive or negative exponential law. If p < 1, then we can 
leave aside the exponential part inside the infinite product for a canonical representation and c = 0. In any case 
tp(t) is factorized as a convergent product of characteristic functions, and hence it is a characteristic function 
due to Levy's continuity Theorem. 

The two mentioned examples differ in terms of the density function. Notice that the density deduced from 
F\(x) can be considered an infinite linear combination of exponential densities, while f2{x) seems to be a mix 
of Laplace densities. The difference between both cases lies in the entire function g{z). On the one hand, gi(z) 
has order 1/2 and all its zeros are simple and positive; as opposed to the second example, where g2 (z) has order 
1 and its zeros are simple but symmetric with respect to the origin. 

The aim of this paper is to give some general results suggested by those and other similar examples. It 
will be shown that when g(z) and h(z) are both entire functions of order p, p' £ (0, 2) respectively, satisfying 
g(0) = h(0) = 1 and with a further condition over their zeros then tp(t) = h(it)/g(it) is a characteristic 
function. We prove that when p, p' € (0, 1), the zeros of h{z) and g{z) are simple and positive, and some 
other additional hypotheses, then the probability measure corresponding to ip{t) has a density f(x) that can be 
written as a sum of exponential type densities, 



This series is a generalized Dirichlet series (see Mandelbrojt lfT6l ) which has very good properties; in particular, 
we prove that it is uniformly convergent on every compact subset of (0, oo), and the cumulative distribution 
function also turns out to be a generalized Dirichlet series. 

When the zeros are simple and symmetric, positive and negative, the existence of a density can be also 
proved for p = 1 and h{z) = 1. Then the density can be written as a series of Laplace type densities. This case 
is important because it covers some elements of the double Wiener chaos (see Janson flU) as the double Ito- 
Wiener integrals where the kernel has symmetric zeros with multiplicity 2; in particular, the Levy area belongs 
to this class. 

We study some examples. In addition to the square of a Kolmogorov law and the Levy area, we consider 
the law of the first hitting time of a Bessel process, whose characteristic function is expressed as a quotient of 
Bessel functions (Kent [12] and Borodin and Salminen O). We also show how this technique can be used to 
invert some Laplace transforms. In our last example we study a particular case of the Heston model used in 
mathematical finance and we prove that the general theory developed in the first part can be applied to it; such 
a study was the starting point of this paper. 

To summarize, the paper is twofold. On the one hand, it shows a way to construct a rich family of charac- 
teristic functions. On the other, given a characteristic function that can be identified as belonging to that family, 
our results provide a procedure to invert that characteristic function. 

2 Construction of characteristic functions 

In this section we identify a rich family of characteristic functions which can be constructed by quotients of 
entire functions. 




1 
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From now on we will consider that a strictly increasing sequence of positive numbers, {a n n > 1}, is a 
sequence satisfying < a\ < a2 < • • • and such that lirrin^oo a n = oo. 

Proposition 2.1. Consider two entire functions g{z) and h(z) of order lying in (0, 1) such that g(0) = h(0) = 
1, with simple positive zeros given by the strictly increasing sequences of positive numbers {a n , n > 1} and 
{b n , n > 1} respectively. Assume that a n < b n for all n. Then tp(t) = h(it) / g(it) is a characteristic function 
of a probability measure on [0, oo). 

Proof. By Hadamard's factorization Theorem (see, for example, Levin lfT3l p. 26]), 

OO / \ OO / 

5 (z) = n(l-f and h(z) = H(l 



n=l v ' n=l 

We conclude from the convergence of both products that 



b„ 



n=l n=l 

Since < a n /b n < 1 it follows that each factor of the above product is the characteristic function of the 
probability measure 

£ <50+ ( 1 "^) £Manh 

where <5o is a Dirac measure at and £xp(a n ) is an exponential law with parameter a n . The result is a 
consequence of Levy's continuity Theorem. ■ 

Remarks 2.2. 

1) The condition g(0) = h(0) = 1 is merely a way to ease the notation, in fact the same result would be true 
if we let g(0) = h(0) / 0. 

2) By means of standard manipulations, Proposition 12. 1 I holds true for h = 1. In this case each factor of (p(t) 
is the characteristic function of the probability law £xp(a n ). 

3) A similar result to Proposition 12.11 is proved in Corollary 9.16 by Schilling et al. lfT9l . where it is also 
deduced that the corresponding probability measure belongs to the Bondesson class, and if h = 1, then 
the probability measure is in the class of convolutions of exponential densities (for the definition of these 
classes see Schilling et al. ITT91 Pages 80 and 87]). 

4) It is easy to show that the restriction of simple zeros in Proposition 12.11 can be relaxed if we let a n and b n 
have the same multiplicity for all n. 

Proposition 2.3. Consider two even entire functions g(z) and h(z) of order lying in (0, 2) such that g(0) = 
h(0) = 1. Let {a n , n > 1} and {b n , n > 1} be two strictly increasing sequences of positive numbers; and let 
{±a n , n > 1} and {±6 n , n > 1} be the simple zeros of g(z) and h(z) respectively; assume that a n < b n for 
all n. Then <p(t) = h(it) / g(it) is a characteristic function of a probability measure on R. 

Proof. By Hadamard's factorization Theorem 

°°/2\ 00 /2\ 

71=1 ^ IL/ n =l v u/ 

due to the symmetry of the entire functions. The rest of the proof follows the same argument as the derivation of 
Proposition 12. II In this setting each factor of tp(t) = h(it)/g(it) is the characteristic function of the probability 
measure 2 2 

^ <5o + - pM £aplace(a n ) , 
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where £aplace(a n ) is a Laplace law with parameter a n . ■ 

As pointed out in Remark I2T21 there is a straightforward generalization of Proposition 12. 3 1 when h = 1. 

The statements of Propositions 12. 1 1 and 12.31 are related to the main results of the following sections. The 
reader will notice that it is possible to state more general results using the same ideas and hence with similar 
derivations. We may state Proposition 12.31 without requesting the functions to be symmetric but with their zeros 
given by the sequences {a n , n > 1} and {b n , n > 1} where a n and b n have the same sign, and \a n \ < \b n \. In 
such a case 

V {t) = = J* ft + (l - - -Y 1 ) e<^) 

9W „^ \ bn V b nJ V <W J 

(assuming d € R) and each factor is the characteristic function of a translation of the probability measure 

For example, we can set ip(t) = Ai(uit) / Ai(vit), where < u < v < oo and Ai(z) is Airy's function. The 
Airy function is an entire function of order 3/2 and its zeros are real and negative (see Katori and Tanemura 
ifTTTO . In this case, although we can state that ip(t) is a characteristic function we will not be able to write its 
distribution function explicitly. 



2.1 Finite convolution of exponential and Laplace densities 

The next two lemmas will prove useful for determining the density of the characteristic functions in Proposition 
I2.1l and l2~3l we will consider a finite product approximation of tp(t) and determine its density so we can take 
the limit afterwards. This section and the one after will develop the first steps of the procedure. 

Let us recall some standard notations. Given two probability measures on (R, £>(R)), Pi and P2, we denote 
by Pi * P 2 its convolution 

Pi*P 2 (B) := [ P 1 (B-y)P 2 (dy) , 

where B £ B(M) and B — y = {x — y,x<E B}. The characteristic function of Pi * P 2 is the product of the 
characteristic functions of Pi and P 2 . Moreover, if Pi and P 2 are absolutely continuous with density fi and / 2 
respectively, then the density of Pi * P 2 is given by the convolution of f\ and / 2 

/+00 
fi(y)f 2 (x-y)dy. 
-00 

The convolution of Pi * • • • * P n (resp. /1 * ■ • • * f n ) is denoted by *j =1 Pj (resp. *y =1 /j). 

Lemma [2~4l is well known in the literature (see, for example, problem 12 of chapter 1 of Feller lT7lD. 

Lemma 2.4. Fix < Ai < A 2 < . . . < X n and define the couple A{n) := Y\2=i an ^ B(k,n) := 
YYi=i (Afc — Ai). Then *? = i£ xp(A,) has density given by 

n —X-x 

/n (x) = (-ir+^(n)^^— - x>0. 

~^ B(i,n) 



The proof of the following result uses an interesting property given by Bondesson (see O) in the context 
of generalized gamma convolutions that, in our setup, states that if Y is a non-negative random variable with 
moment generating function My(u), for it in a neighborhood of 0, and T is a centered normal random variable 
with variance 2, independent of Y, then the random variable X := VY T has moment generating function 
Mx{u) = My(u 2 ). The proof is completed in the line 



M x (u) 



E[e u ^ T ] 



E 



E[e uVYT /Y] \ =E[e u 



M Y (u 2 ) . 



4 



In fact, the next result is a consequence of the previous lemma and will be useful for deriving the density of the 
characteristic function of Proposition 12.31 

Lemma 2.5. Fix < Ai < A2 < . . . < A n and define E(k, n) = YYi=i (^1 ~ )• Then *™ =1 £aplace(A,) has 
density given by 

(_y\n+l " P -Ai|x| 

/ „ w = L^^ ( „ ) g__ xeR . 

Proof. Let us consider the characteristic function 

which corresponds to a random variable, Y, that is the sum of n independent exponential random variables with 
parameters {A|, 1 < j < n}. By the previous lemma, the density of Y is 

w) = (-irA.)Eip »^»- 

Let T be a centered normal random variable with variance 2, hence the random variable X := VYT has 
characteristic function 

which corresponds to the sum of n independent Laplace random variables with parameters {Xj, 1 < j < n}. In 
order to find the density of X, consider the pair (X, T) and compute its marginal density by means of a change 
of variables to the pair (Y,T). The Jacobian determinant is -tU-*-, which means that the change of variables is 
an almost everywhere diffeomorphism of the plane. Fix then 

/•ca rO 

fx(x) = l x > f\ x ,T)(x,t)dt + l x<0 f(x,T)(x,t)dt 

JO J-00 




as required. ■ 

2.2 Finite density approximation of the characteristic function 

We derive here the density function for a finite approximation of the characteristic functions of Proposition [2J] 
andO 
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Lemma 2.6. Let < a\ < (12 < ■ ■ ■ < a n and < 61 < 62 < • • • < b n such that aj < hi for 1 <i <n. Write 
either 

(a) g n {z) = n" = i (1 - z/cti) and h n (z) = Y\™=i (1 ~ z /h), or 

(b) 9n(z) = m=i (1 - *7a?) and h n (z) = EKLi (l - z 2 / 6 ?)- 

Let ip n (t) be the characteristic function <p n (t) = h n {it) / g n (it). Then the corresponding law ofcp n (t) is 



n 

\i=i 



04 



5q + > 



w/zere /x„ jj a finite measure on (0, 00) - case (|3p - or on M \ {0} - case © - w/f/z density given by 



(tap 



dx 



x > or 



h n {aj) - aj \ x \ 



x^O 



Proof. We will prove the result for the case ©; case © is similar. Since the convolution of two finite measures 
is absolutely continuous when one of them is, we abuse of the notation and write 



1 



£xp(aj) ] (x 



to denote the probability density function of the measure which is a mixture of a delta measure and an expo- 
nential distribution. Let us denote by I l n the set of all subsets of {1, 2, . . . , n} of cardinal 1 < i < n. Then 



fir Uw+E £ 



\i=l 
f n 



i=1 J&I™-* 



n 

fceJ 
a " x a k 



bk 



n 1 



fceJ c 

n 



OA- 



£ 



~ a k x 



vt=l V fc=l llr=l I 1- ^ I i=l j 6/ «- 



u=l 



r=£k 



a k x 



fceJ c TlreJ c \{fc} 



1 _ Ofc 



reJ 



fc=l 



r^k 



1 a fc 



n 1 



r=l 



+ 



r£J c 



n £)*<«> 



u=l 



n n n ( 1 — 5& 

^ in^i(i-t 



where we have used Lemma [2~4l ■ 

Remark 2.7. If h = 1 in Lemma I2T61 then the law of <^ n (i) is just fi n . 

Lemma 12.61 shows that the distribution has an atom at the origin. This causes a handicap towards apply- 
ing the conventional limit theorems to the expressions of the densities therein. Section 14.21 will give further 
assumptions and criteria to overcome this difficulty. Finally, notice that the key point to this end is determining 
the behaviour of 

Iim — - — . 

n— j-oo dx 
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3 Fundamental lemmas 



The following lemmas are essential for the purposes of the paper and they give the technical results to allow us 
to use the standard convergence results on the expressions of Lemma [Z6l Before that, let us remark on a key 
property of an entire function: the order of an entire function is greater or equal to the exponent of convergence 
of its zeros (see Titchmarsh EOl p. 251]). That means, let / be an entire function of order p and {a n , n > 1} 
its zeros, then 



p > inf{a > : V] \a n \ a < oo} . 



n=l 

In particular, Yl^Li | a n | < oo for /3 > p. The equality holds if / has a representation as a canonical product, 
in fact Proposition 12. 1 1 and 12.31 give such representation for h and g. 

Lemma 3.1. Let g{z) be an entire function of order p G (0, 1) such that g(0) = 1, with simple positive zeros 
given by the strictly increasing sequence of positive numbers {a n , n > 1}. Let h(z) be another function that 
either 

(a) h(z) = 1 or 

(b) h(z) is an entire function of order p' € (0, 1) such that h(a n ) ^ Ofor all n. 
Then, for every x > 0, the series Y2n>l g^(a \ e ~ a " x is convergent. 

Proof. We will prove the lemma for case (b); case (a) is similar and indeed easier. Consider the closed contour 
D{R) in Figure[TJ where R > 0, 9q G (0, 7t/2) and R ^ a,j for all j. Fix x > 0. The function e~ xz is 



/ 

/ 

/ \ 9 o 



D'(R) / 
/ 

/ 

/ 

/ 



/ 



s 

s 



\ 



\ 

\ 

\ 

D 2 (R) S> 



\ D 3 (R) 

\ 

\ 

\ 
\ 
\ 
I 



j n + 1 




Figure 1. Dashed line is the contour 
D(R) of integration. 



Figure 2. The function \g(z)\ is 
bounded below in the region ex- 
cluded from the gray circles. 



analytic in the region bounded by D(R), except at the zeros a,j that lie in the interior of D(R). The Residue 
Theorem allows us to write 



1 

2vri 



d{r) g(z) 



dz 



E 



h(a 



g'(aj) 



J / g— ajX 



Our objective is to prove that there is an increasing sequence, {R n , n > 1}, towards infinity such that the limit 



Hz) 



h ™ T D{Rn) g{z) 



dz 



(2) 
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exists and is finite. The key point of the proof is that an entire function of order p < 1 can be bounded below in 
the region excluded from small circles around its zeros. Since the zeros of g(z) are located on the positive real 
axis, convenient bounds of 1/|^(^)| can be obtained. Specifically, the proof is based on: 

1. For every e > 0, g(z) = 0(e e ' 2 ') and h(z) = ©(e 6 ' 2 '). This is due to Yl n l a nl _1 < °°> an ^ hence we can 
apply Titchmarsh's issue 15 E01 p. 286]. 

2. As a consequence of the previous statement, for every e > there is an increasing sequence {R n ,n > 1} 
such that lim ra R n = oo and \g(R n e %9 )\ > exp{— eR n } uniformly on 9 £ [0, 2tt]. See Titchmarsh l20l p. 
276, it. 8.75]. 

3. For every e > 0, there is 9q 6 (0, vr/2) and vq > such that 

\g(re ±ie °)\ > exp{— er} , for r > tq. 

This is deduced from Titchmarsh GUI p. 273, it. 8.71]. Take straight lines through the origin with angles 9q 
and —9o that do not intersect with any of the circles with centers aj and radius 1/oj (for aj > 1), see Figure 
[2] Titchmarsh proves that for all e' > exists r' such that if r > r' we have \g(z)\ > exp{— r p+e } in the 
region excluded from these discs. Take e' satisfying p + e' < 1 and thus for sufficiently large r we have 
r p+e £r Choose ro a slightly larger than r' in order to obtain the claimed inequality. 

Now we are ready to complete the proof. Fix e > such that 2e < x cos 9q, where 9q is the angle depending 
on e such that property [3]is fulfilled. Denote by D 1 (R n ) and D 2 (R n ) the lower and upper straight segment of 
D(R n ) while D 3 (R n ) stands for the arch. All paths are considered with the corresponding orientation. Then 
we divide the integral into three parts : 



D 



(Rn) 9{z) 



dz 



I +/ + / 

JD 1 (R r ,) JD 2 (Rr,) JD- 



'Di(^n) JD^(R n ) JD3(R n ) 

We first consider the integral over D 3 (R n ). Due to points [J and |2l we can bound the module of this integral 
in the following way: 



f Mf) 



dz 



<Rn 



h{R n e ie ) 



g{R n e ie ) 

< ft e ~xR n cos 6»o 



-xR n COS t 



de 



rd 


h(R n e i0 ) 


l-0o 


g{R n e ie ) 



de < KeoR n e- R ^ xcose °- 2e \ 



for R n large enough. This goes to zero as n — > oo. 
The integral over D 1 (R n ) can be parametrized as 



D 1 (R n ) 



We claim that 



dz 



h(re 



-i6 



R„ 



hire 



-iQo s 



g(re~™°) 



dr . 



-ie 



dr < oo. 



g(re~ id o) 

According to observations [1] and [3j there is ro > depending on e such that 



(3) 



h(re~ ieo ) 



dr 



ro 



hire 



-^6»o^ 



g{re 



-xr cos Qq 



Since gire 
hence 



g(r e -i0o) 

~ ld °) / for r > 0, it turns out that the function 



dr < K 



-r(xcos6» -2£) g r 



< OO. 



hire 



-iOo^- 



gire~ ie °) 



-,—xr cos 9o 



is continuous on [0, ro] and 



dr < oo. 



o \g(re- ie °) 

Adding up the two upper bounds we obtain (f3]>. For D 2 (i? ra ) the computations are equivalent and the limit © 
exists and is finite. ■ 
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Lemma 3.2. Under assumptions of Lemma li.il the series 



converges absolutely for x > 0, ant/ ?/ze convergence is uniform on any compact subset 0/(0, oo). 

Proof. The above expression is a generalized Dirichlet series. As a consequence, if (01) is convergent for some 
xq, then it is so for all x > xq and the convergence is uniform on every compact subset of the half-line (see 
Mandelbrojt [16, p. 9]). Therefore the second part of the lemma follows from Lemma [3TT1 Each general 
Dirichlet series has associated an abscissa a c of convergence and an abscissa a a of absolute convergence. In 
general these two values are not equal but their distance is bounded (see Mandelbrojt [16, p. 1 1]) by 

Inn 

< a a — o~ c < lim sup . 

From the previous observation we know that a c < 0; we claim that the above limit is zero to derive the first 
part of the result. The claim is proved through Jensen's formula which gives the relationship n(r) = 0(r p+e ) 
for the entire function g(z), where e > and n(r) stands for the number of zeros with norm less than or equal 
to r (see Titchamrsh |[20l p. 249]). Choose e such that p + e < 1, then n < Kan +£ for some positive constant 
K. Finally 

Inn < n- 1 < K< +£ - 1 , 

and hence a c = a a . ■ 

Remark 3.3. Straightforward manipulations lead to generalizing Lemmas 13.11 and 13.21 for the case where the 
simple zeros of g(z) are {±a n , n > 1}, but still of an order lying in (0, 1). 

The next result somehow extends case (a) of Lemma [3TT1 but with a penalization of extra hypotheses. As 
will be shown in the examples, this particular case is also of great interest. 

Lemma 3.4. Consider an even entire function g(z) of order 1 such that g(0) = landg(z) = 0(e A \ z \) for some 
positive constant A. Let {a n , n > 1} be a strictly increasing sequence of positive numbers; let {=ba n , n > 1} 
be the simple zeros ofg(z); assume that the following limit exists: 

Tl 

lim — = 5 > . (5) 

Then, for every x > 0, the series Yl n >i g '( a ) e~ a " x i^ absolutely convergent and the convergence is uniform 
on any compact subset of(0, oo). 

Proof. It is clear from (f5]) that a a = a c , thus if we prove the plain convergence of the series in (0, oo) the result 
will follow. The proof is very similar to that of Lemma l3.1l but the fact that the order of g{z) is 1 demands some 
modifications. We will consider as before the closed contour D(R) in Figure [TJ where R > 0, Q$ G (0, 7r/2) 
and R ^ a,j for all j. Our objective is to prove that there is an increasing sequence, {R n , n > 1}, towards 
infinity such that 

lim / -!- e~ xz (6) 

n J D(R n ) 9(Z) 

exists and is finite for a fixed x > 0. Split the above integral as in Lemma [37T1 into D(R n ) = D 1 (R n ) U 
D 2 (R n )UD 3 (R n ). 

The function g(z) admits an expression as a canonical product given by 



»(*) = n 



n=l 



Z 2 



Or 
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Define the function 



OO , s 

^) ==n( i - i) 

n=l V 



of order 1/2 and notice g(z 2 ) = g{z). Fortunately, the function g(z) inherits the good properties of g(z). To 
start with, since g(z) has order 1/2, for every e > there is a sequence R n /* oo such that for all G [0, 27r] 

where A%(r) = max|| 2 | =r } |g(z)| (see Titchmarsh ll20l p. 275, it. 8.74]). Since g(z) is even, we deduce that 

— 1/2 

there is a sequence R n := i?„ oo such that for all 9 G [0, 27r] 

| 5 (i?„e ie )| > M s (i4) _e • 
Given that #(z) = C(e A l z l), for i? n large enough, 

1 



g{R n e ie ) 

Now consider e such that Ae < x cos #o> an d the bound 

1 



< Ce eAHn uniformly on 9 G [0, 2vr]. 



e~ xz dz 



<Rr, 



-e 



-xR n cos ( 



de 



< ^ g— xRn cos 9 



[00 


l 


Le 


g(R n e ie ) 



d9 < K6 R n e- Rn{xcos9 °- A£ \ 



for R n large enough. This goes to zero as n — > oo. 

In order to bound the integral over D 1 (R n ) consider the following Remark in Levin lfT3l p. 82, eqn. (2')]: 



ln\g(re id )\ = ? 1/2 ir5 sin(0/2) + 



Q(rV 2 ) 
trin(0/2) 



for 9 G (0, 27r) and as r — >• oo; that asymptotic result is based on the fact that g is a canonical product of order 

1/2. Notice that g(re ld ) = ^(r 1 ' 2 ^' 2 ) := g(re ld ) = g(re i( - 9+7T ^), hence the above asymptotic equation is 
translated into g(z) as 



In \g{re M )\ = rir5\ sin(0)| + 



o(r) 
sin(0)| 



for (9 G (0, 2?r) \ {vr}. Therefore 
1 

/Di(R„) 9(2) 



e~ xz dz 



< 







-ir cos 9q 



-r(a cos e +7nS| sin(0 o )| + r | il^yj ) rfr 



where ~ means that both converge or diverge together. Clearly, the last integral is finite. The same derivation 
is valid for the integral over D 2 (R n ) and © holds. ■ 



4 Construction of density functions 

This section will derive the density function associated to the characteristic functions of Proposition 12.11 and 
Proposition [23] under the restrictions of Lemma [3^41 

Lemma 12.61 shows the density function for a finite product approximation of the characteristic function. 
From Levy's continuity theorem, this means a convergence in distribution of the associated probability mea- 
sures. We first show a pointwise convergence of the distribution function to finally obtain the convergence of 
densities. In order to avoid Dirac's delta measure at 0, instead of using the distribution function F(x) we work 
withF(x) := 1-F(x). 
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Theorem 4.1. Under assumptions of Proposition |Z71 the probability measure on [0, oo) corresponding to the 
characteristic function ip(t) = h(it)/g(it) is absolutely continuous on (0, oo) with (perhaps defective) density 
given by 

/(*)=-f;4H e ~ a ^ x>o - 

Proof. Recall h n (z), g n (z) and f n (t) from Lemma 12761 case fla]). Denote by F n (x) the distribution function 
corresponding to the characteristic function (p n (t). It is clear that ip n (t) — > ip(t) pointwise, where <p{t) is 
defined in (fl]). Since a n / for all n, it follows that ip(t) is continuous at 0. Therefore there exists a distribution 
function F(x), such that 

lim F n (x) = Fix) 

n— >oo 

for all x where F(x) is continuous. From the expression of the density of F n (x) we deduce that 

F n {x) = -j2±4^e~^, *>0. 

^ a i9n{ a i) 

Our objective is to prove that 

lim£^^e-- = f;4^-- x>0. (7) 

The proof of (0 is done using the dominated convergence theorem. First notice that for alii > 1, 

lim g' n (ai) = g'{ai). (8) 

n— >oo 

This is proved in the following way: On the one hand, 

lim g' n (ai) = -- lim II f 1 " -) = - 1 II f 1 " -)• (9) 



On the other hand, for z ^ aj , 



a, n->oo -*- v a, / aj 

i=i J j=i 



1 - z/aj 

The function on the right hand side is entire, so by analytic continuation, 



T^-n (»--)■ < io > 



Further, the function l/g(z) has a simple pole at Oj. Hence, 



lim (z - aj)— !— = Res(l/q,aj) = 1 (11) 



Combining ©-(QT|> we obtain ®. 

Fix x > 0, set A = {aj, i > 1} and denote the counting measure by V. By ©, 

ra„(y) 2/0 (y) 



n 



Now observe that for n > j 



K+i{aj) _ h n (a,j) l-aj/b n 



+1 



5n+i(°i) 

and since a n+ \ < b n+ i we conclude that 



9n(. a j) 1 - aj/a n+1 ' 



hniflj) 


< 


h n +i{cLj) 


< 


h(dj) 


g'rMj) 




9n+l(«i) 




9'(aj) 



(12) 



Finally, let <I>(y) be the function 

Kv) 



$(y) := 



Then 



K{y) - 



!)■>' 



ya'niy) 

By Lemma I3T21 



< 



yg'(y) 
K{y) 



L [o,i 



(y)e 



yy'niy) 



l[o,i) (l/K 



-yx 



+ 



+ 



Kv) 



g'(y) 
h n (y) 



1 [i,oo)(2/)e" 



l$|dr 



i: cii<l 



h{ai 



a%sf{<H) 



9n(y) 

aiX _|_ ^ y 



i,oo)(y)e 
h(di) 



-yx 



-yx 



y>0 



< Hy), 



y>0 



i: a;>l 



g'(ai) 



Hence we can apply the dominated convergence theorem to prove ©. It follows 

F(s) = \ r .r>0 



\ - h(aj) 
~ a»5'(oi) 



Also from Lemma [X2l we check that the set of continuity of -F(x) is (0,oo). Moreover, from the uniform 
convergence on compact sets of (01), 



^7? 



x > 



as required. ■ 

Remark 4.2. Standard manipulations show that Theorem 14. H is also true for h = 1. 

Remark 4.3. Notice that inequality (fl2l also holds true for Lemma 1231 case ©. Thus, Theorem 14. II for such 
g{z) and /i(z) with orders lying in (0, 1) still holds and concludes that the limiting density is 



/(*) 



/i(a n ) 



^ — \ n 



-a n \x\ 



Moreover, the same proof of Theorem 14.11 will be valid for functions g(z) and h{z) of order lying in [1, 2), 
provided that a a < for the series (0]). 

The next result is an extension of the preceding theorem on the existence of a density function. One would 
like to generalize such a result for characteristic functions ip(t) = h[it)/g(it) where g(z) and h(z) have zeros 
{±o n , n > 1} and {±6 n , n > 1} and orders lying in [1,2). Entire functions of order greater than or equal 
to one are much more difficult to treat than those that have order less than one. Therefore we will restrict the 
study to the setup of Lemma [3~4l As pointed out in the above remark, we only need the absolute convergence 
of the general Dirichlet series to follow the same proof of Theorem 14. 1 1 and obtain the density 



/(*) 



-a n \x\ 



X^O. 



(13) 



However, we would like to present a different proof, Theorem 14. 4[ which shows that both problems, the expo- 
nential and Laplace convolutions, are the head and tail of the same coin. It is worth remarking that the random 
variables corresponding to that case are in the homogeneous second Wiener chaos (see Janson |8, chap. 6]), 
so our result is a step forward in the study of the densities of such an interesting space of random variables; in 
particular, this case includes the Levy area. 
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Theorem 4.4. Under assumptions of Lemma \3.4\ the probability measure on K corresponding to the charac- 
teristic function ip{t) = 1/ g(it) is absolutely continuous on R\ {0} with (perhaps defective) density given by 

(123. 

Proof. Recall the expressions of g(z) and g(z) defined in the proof of Lemma l3.4l Notice that 

2a0{a 2 j ) = g'{a 3 ) . (14) 

Due to Proposition 12. 1 1 and Theorem 14. II ip(t) = l/g(it) is a characteristic function of a non-negative random 
variable, denoted by Y, with density 

Fix x 7^ and proceed heuristically as in Lemma [231 using the afore mentioned Bondesson argument and (fl4l 
to obtain the required expression (fl"3"l) . To make the argument accurate we need to apply Fubini's theorem since 
jV(x) is an infinite series. It turns out that the absolute convergence of (0]) allows us to use Fubini's result. 



4.1 Infinite convolution of exponential and Laplace densities 



For the sake of completeness we will rewrite Remark [472] and Theorem 14 .4 1 in a way that will extend Lemmas 
2 .4 1 and 1231 Moreover, we prove that in this setup the resulting density is continuous in [0, oo) or IR. 

Following Wintner, for a sequence of densities, {f n , n > 1}, we will say that *^ = if n is a convergent 
infinite convolution if the product 



where 



n=l 



Atx 



fn(x)dx , 



is uniformly convergent in every fixed finite i-interval. 

Proposition 4.5. Let {X n n > 1} be a strictly increasing sequence of positive numbers such that for some 
p £ (0, 1), Yl n >l < °°- Then ★^L 1 fxp(A n ) converges to a continuous density on [0,oo) which can be 
written as 



K n=l 



£xp(A n )(x) = lim y~] 

n— >oo — ' 



i=l 



B(i,n) 



8=1 



k=l 



-i 



e~ XlX , x > 0, (15) 



and 



*~ 1 5xp(A n )(0) = 0. 



Proof. We will first point out that the infinite convolution is convergent. To this end we will use a very useful 
result from Wintner [21] that ensures the convergence of the infinite convolution if 



oo „ c 

^M n < oo where M n = E[|£xp(A n )|] = / 
n=l 



\x\\ n e An2: l[ 0iOO )(ix . 



Clearly M n = A^ 1 and the condition to guarantee the convergence of the infinite convolution is fulfilled. 
Moreover, if one term of the infinite convolution has a continuous density of bounded variation, then so has the 
infinite convolution; and the continuous function *^ =1 £xp(A„) tends, as m — > oo, to the infinite convolution 



13 



uniformly in every bounded range (see Wintner |[22l "). It suffices to notice that £xp(Ai)*£xp(A2) is continuous 
and of bounded variation. Hence, for all x > 0, 

*™ =1 SMK)(x)= lim E ^R^ll^ ^ 

1=1 y ' ; 

To prove the second equality in (TTSb . let 

00 / \ 

^=n(i-^)- 

Then g(z) is an entire function of order less than 1 and apply Theorem 14. 1 1 to x > 0, and the computations 
done in the first part of the proof of that theorem. The value of the density at zero is guarantied by Wintner |[22ll 
result about the continuity on R of the density. ■ 



Note that the series in the right hand side of (fT5l) may be divergent at x = 0. 

Proposition 4.6. Let {X n n > 1} be a strictly increasing sequence of positive numbers such that Yl n >i ^n 2 < 
00. Then *^ 1 £aplace(A n ) converges to a continuous density on R which can be written as 

C=i^aplace(A ri )(x) = lim V ^ '—. — — — e~ Ai ' :c ' , x G R . 

n-s>oo^^ E(i,n) 
i=l v ' ; 

Proof. As in Proposition [43J we start by showing that the infinite convolution is convergent, lessen and Wintner 
show in J9l that if 

00 00 
Y,M l n and Y, M n 

n=l n=l 

are convergent, where M\ and M% are the first and second moment of £aplace(A n ) respectively, then so 
is the infinite convolution. One can check that M\ = and M% = A^ 2 to obtain the convergence of the 
infinite convolution. Therefore we can apply Wintner |[22l to obtain the proposition since a Laplace density is 
continuous and of finite variation. ■ 



4.2 Existence of densities 

In this section we give three different conditions that guarantee that the probability measure corresponding to 
h(it)/g(it) is absolutely continuous. 

Proposition 4.7. Under the hypotheses of Theorem \4.1\ ifh(z) = 1, then the probability measure correspond- 
ing to l/g(it) is absolutely continuous and its density is continuous. 

This result follows from Proposition [43] 

In some cases it is known that the probability measure corresponding to h(it)/g(it) has no atom at zero. 
Since that probability measure is concentrated on [0, 00), this suffices for the existence of a density. 

Proposition 4.8. Assume the hypothesis ofTheorem \4. 1 \ and denote by fi the probability measure corresponding 
to h(it)/g(it). If /u({0}) = 0, then \i is absolutely continuous. 

The following lemma gives a sufficient condition in order to apply a classical criterion for the existence of 
a continuous density. 
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Lemma 4.9. Let g(z) and h(z) be two entire functions of order p £ (0, 1), both with non-zero positive simple 
real zeros, and the zeros of g(z) different from the zeros of h(z). Denote by n g (r) (respectively n/ l (r)J the 
number of the zeros of g(z) with module less than r. Assume the existence of the limits 



n„(r) 
lim > 

? — >oo tP 



and 



5' 



lim > 0, 

r— >oo tP 



with 5' < 5. Then 



/•oo 


h(it) 


1 — oo 


giit) 



dt < oo 



Proof. Levin ( lfT3l p. 82, eqn. (2')]) proves that if 9 € (0, 2tt) and r — > oo then 

ir5r p cos (p(6 — ir)) o(r p ) 



log| 5 (re ie )| 

Sin(7T/9j 

The analogous holds true for h(re ) with 5' instead of 5. 
Split integral (TT6l) into two parts as 



+ 



h(it) 



giit) 



dt 



g(te i7T / 2 ) 



dt + 



sin(0/2) 



h(t^l 2 ) 



(16) 



(17) 



g(te i37r / 2 ) 



dt. 



(18) 



Due to (TTTT ) there is vq > and C > such that for r > tq we get 

/i(te^/ 2 ) 



g(te in / 2 ) 



Hence 



exp{— cos(/97t/2)7t(5 — 5') csc(7rp)r p + o(r p )} 
exp{-r p (cos(p7r/2)7r(5 - 5')csc(i:p) + o(r p )/r p )} < e 
h(te™l 2 ) 



CrP 



dt < 



-CrP 



dr < oo 



g(te i7T / 2 ) 

since the last expression can be reduced to a convergent gamma integral. For the integral over [0, ro] there is 
a straightforward bound using the same sort of derivations used in Lemma |3?T1 The other integral on the right 
hand side of (fT8l ) is bounded in a similar way. ■ 
With the notations of this lemma, 

Proposition 4.10. Under the hypothesis of Theorem \4.1\ and assume the existence of the limits 



n a (r) 
lim > 



and 



6' 



lim — ^ > 0, 



! — >oo tP r— >oc tP 
with 5' < 5. Then the probability measure corresponding to h(it) / 'g(it) has a continuous density. 



5 Examples 

We will now see different situations where we can apply the results obtained in the previous sections. Most 
results are known, but here we get them all using the same technique. 



5.1 Levy area 

Let (p(t) = sech(£T) be the characteristic function for the Levy area (see Levy lfT4lO . where the time component 
of the process varies in [0, T]. Here we set h = 1 and g(z) = cos(zT), where g(z) has order 1. It is clear that 
{±(2/c — l)7r/2T, k > 1} and {±(— l) fc /T, k > 1} are the poles and the residues of l/g(z) respectively, thus 
tp(t) fulfills Proposition 12.31 Moreover, g(z) is even and of exponential type, thus we apply Theorem 14.41 to 
obtain 

k=l 

Proposition 14. 61 ensures that the density function is continuous in R. 
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5.2 The first hitting time of a Bessel process 



The second example is a characteristic function obtained from 

k(z) = z~ vl2 2 v l 2 T(v + i)J„(V2z) , 

where J u (z) is the Bessel function of first kind and order v > — 1. Consider < u < v < oo and let 
h(z) = k(u 2 z) and g(z) = k{v 2 z). The probability measure of the corresponding characteristic function 
(p(t) = h(it)/g(it) describes the first hitting time of the point v by a Bessel process of order v that starts at u, 
see Kent Ifl2l . and can be expressed as 



;v\v J u {uy/2it) 

fit) - 



^uJ J u (vV2it) 
From the Taylor expansion of 

I I z 

v 7 n>0 V ; 

we can deduce the order of k(z), since the order of an entire function is 

nln(n) 



p = limsup 



ln(l/|c n | 



where c n are the coefficients of the Taylor expansion, see Levin [13, p. 6]. Due to the Stirling formula the 
above limit for h(z) and g{z) is 1/2. Denote by {j v j., k > l}the positive zeros in order of magnitude of the 
Bessel function J v [z) and by = j 2 k /(2v 2 ) the zeros of g{z). Finally we can apply Theorem 14. II to obtain 
the density 

/(X) = g U»J v+1 Uu, k ) ^ X>0 ' (19) 



where we have used 



^-Jy{z) = -J V+1 (Z) + -J V (Z) 

az z 



Equation ( fl9l ) is also derived by Borodin and Salminen (3l p. 387]. Notice that the distribution function is 
absolutely continuous in [0, oo) since the probability distribution gives no mass to {0} due to the continuous 
paths of the Bessel process. 



5.2.1 Exit time from a n-dimensional sphere by a Brownian motion 

Let T n denote the random variable of the total time spent by an n-dimensional Brownian motion starting at 
inside the sphere 5 n_1 (r) of radius r > and n > 3. Let P n denote the first exit time for a n-dimensional 
Brownian motion starting at from the sphere S" 1 ~ 1 (r) for n > 1. Ciesielski and Taylor [4] show the remark- 
able equality of the distribution functions of T n and P n -2 for n > 3, and they derive the distribution function 
of T n using methods developed by Kac [TOl . They first compute the solution for n = 3 and then make a guess 
for the general framework. Finally they compare the result with the distribution of P n which was computed by 
Levy. 

We can use Theorem 14. II to derive the density function of T n since Ciesielski and Taylor give its character- 
istic function. In fact they establish the following result 

\v-l 



(rV2zy- L fr / 2r 2 z 



E[e zTn ] = =■ = TT 



r(u)J u ^(rV2z) f = \ \ j 2 _ h 
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where v = (n — 2)/2. This is a particular case of the previous example, where we let h(z) = 1 and g(z) = 
k(r 2 z). Use the same arguments as before to derive the density function 

f(x)= V 'r 1 * , e-^ x x>0. 

Notice that T n is absolutely continuous due to Proposition 14.7 1 and we can set /(0) = to obtain a continuous 
density. 



5.2.2 The area under a squared Bessel bridge 

The characteristic function of T n is easily representated when v = n + 1/2 for n G N. Let v = 3/2 and r = 1 
to obtain the expression 

m = ^^21 
991 ' sinh(7^2ri) ' 

Revuz and Yor lPT8l p. 465] give the Laplace transform of the area under a squared Bessel process starting at 
any point and arriving at zero. It turns out that the above characteristic function corresponds to a Bessel process 
of order 2 starting and arriving at zero. Such functions appear recursively in the literature and many authors 
have studied them, for instance [ Q and ifTTll . 

The factorization of the characteristic function is particularly easy, and one can check the identity 

— -frfi-— 

where the residues of l/g(z) are {(— l) k 7r 2 k 2 , k > 1}. Finally, the density function can be written as 

Ej^it- 1 )*^ 2 *^-* 3 * 3 */ 2 for x > 

for x = ' 

as stated in (Tj|. This characteristic function also corresponds (with a change of parameters) to the square of a 
Kolmogorov law. The corresponding disttibution function is 

00 

tf 4 = i + 2^(-l) fc e -^ 2fc2:c/2 forx>0, 
fc=l 

which was obtained by Dugue 0. 





5.3 Inverse Laplace transform 

Theta functions and related expressions have proved useful for manipulations of functionals of Brownian mo- 
tion. For instance, Borodin and Salminen use the inverse Laplace transform of 

, . vsmh(u\f2z) 
usmn^v 2z) 



which turns out to be 



v ^\ v — u + 2kv ■■ ■ " ' -' ■ ■ 

4 MW= ;£,-5s?i5- e * »>»■ < 20 > 

Let us consider (p(t) = h(it)/g(it), where g(z) = sin ^^?^ and similar for h(z). Standard manipulations 
lead to the computation of the sequences {— -^r, k > 1} and {(— l) fc+1 ^f sm (^ 7r ) > k > 1}, which are the 
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poles and the residues of h(z)/g(z) respectively. As pointed out in the previous example, both functions h(z) 
and g(z) are entire functions of order 1/2 and ip(t) satisfies the hypothesis of Proposition 12.11 Notice that the 
poles are negative and hence we have to apply a generalization of Theorem 14.11 for the negative case, this is 
straightforward and we obtain the expression 

/(y) = E(-l) fc+1 - S m(-A:vr)e^ y<0. 
* — ' u \v / 

k=l 

Since we consider the Laplace transform we need to change the sign of y in the above expression and consider it 
in the range (0, oo). After doing that, it turns out that the above series is equal to ([201 due to Poisson summation 
formulae. Moreover, we are able to say that the density function f(y) is continuous in [0, oo) due to Proposition 

14. 101 since n g (r) = , where [x] stands for the integer part of x and similar for n/ l (r). 

5.4 Heston density function 

The authors proved in [5] that the density function of the Heston model is C°° and can be expressed as an 
infinite convolution of Bessel type densities. We give here another expression in a particular case. In fact the 
search for such an expression for the general case was the starting point for the present paper, thus we believe 
it is worth recording it, even though it is only a partial result. 

The Heston model for the log-spot is driven by the following system of stochastic differential equations 

dX t = ~V t dt + y/V t dZt 
dV t = a(b-V t )dt + c v / V t dW(t)_ 

where a, b and c are real positive constants. The processes W and Z are two standard correlated Brownian 
motions such that (Z, W)t = pt for some p 6 [—1, 1]. For the particular case of interest we set 2ab = c 2 and 
consider the volatility process V to start at 0. Then the complex moment generating function of log-spot is 

z(x ~pct/2) a z(x -pct/2) 

E\e zXt ] = - - = - 

L J cosh(P(z)) + 2zg£ S mh(P(z)) g(z) ' 

where xo is the initial point for the process X and P(z) = \J {a — cpz) 2 + c 2 (z — z 2 ). The term e z ( x o-pct/2) 
is merely a decentralisation term, while the function g(z) is an entire function of order 1/2 with all zeros being 
real (see [5]). The function g(z) has negative zeros for p = —1, while for p = 1 and a > c all zeros are 
positive. In any case we can apply Theorem 14.11 to obtain an expression of the density function as a series of 
exponential type densities. Obviously, the zeros of g{z) must be computed numerically. 
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